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ABSTRACT 


Sound  fields  in  wedge-shaped  ocean  layers,  modeling  con¬ 
ditions  on  the  continental  shelf,  have  been  studied  at  the 
Naval  Postgraduate  School  in  the  last  few  years  using  the 
method  of  images.  These  studies  are  carried  further  in  the 
present  work.  The  method  is  implemented  in  different  en¬ 
vironmental  conditions.  This  thesis  examines  the  influence 
of  several  parameters  on  the  sound  field  for  downslope  pro¬ 
pagation  in  a  wedge-shaped  fluid  of  speed  of  sound  C2  over- 
lying  a  slow  bottom  of  speed  of  sound  c^.  On  the  basis  of 
qualitative  and  semi-quantitative  analysis  of  the  behavior 
of  the  pressure-depth  profile  for  various  geometrical  and 
physical  parameters,  we  can  conclude  that: 

1.  A  defined  distance,  the  "characteristic  distance" 
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I.  INTRODUCTION 

A.  SOUND  PROPAGATION  IN  SHALLOW  WATER  CHANNEL 

Experimental  investigations  of  sound  propagation  in 
shallow  water  channels  have  been  done  by  several  investi¬ 
gators.  Shallow  water  propagation  is  of  interest  because  of 
the  applications  to  coastal  defense.  These  investigations 
are  expensive  and  time  consuming.  The  use  of  a  computer 
model  should  provide  a  relatively  inexpensive  alternative  to 
observation . 

One  of  the  techniques  uses  normal  mode  theory.  The 
normal  mode  theory,  which  was  introduced  and  developed  by  C. 
L.  Perkeris  [Ref.  1],  gave  the  exact  solution  in  water  of 
constant  depth.  Further  development  of  normal  mode  theory 
was  made  by  L.  Brekovskikh  [Ref.  2],  who  initiated  pressure 


i 


as  an  integral  involving  Bessel  functions  and  solution  of 
the  normal  mode  equation.  Another  theoretical  approach  to 
sound  propagation  in  a  horizontally  stratified  ocean  of  con¬ 
stant  depth  is  given  by  the  method  of  multiple  scattering 
[Ref.  1].  With  this  method,  all  the  previous  theories  can 
be  simplified  by  conversion  into  an  asymptotic  form  which  is 
valid  when  the  acoustic  wavelength  is  small  compared  to  the 
distance  over  which  the  sound  speed  varies  appreciably. 
These  theories  agree  with  the  laboratory  experiments. 
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For  a  water  channel  with  a  small  bottom  slope,  the  sound 
field  may  be  expressed  approximately  in  terms  of  adiabatic 
normal  modes.  To  facilitate  prediction,  R.D.  Groves,  Anton 
Nagl,  h.  Uberall,  and  G.  L.  Zauer  [Ref.  3]  modeled  a  wedge- 
shaped  isovelocity  ocean  with  a  linearly-sloping,  perfectly- 
rigid  ocean  floor  using  adiabatic  normal  modes.  For  a  pene¬ 
trable  bottom  the  normal  mode  description  fails  when  modes 
propagating  upslope  encounter  the  "critical  depth"  (Hc), 
defined  as  the  depth  where  the  associate  mode  changes  from 
fully  trapped  within  the  water  channel  to  radiating  energy 
into  the  bottom  (cut  off)  [Ref.  4-6].  The  parabolic  equa¬ 
tion  can  be  used  to  explain  the  mechanism  of  sound  energy 
radiation  into  the  bottom  [Ref.  7,8].  Such  an  eauation  was 
studied  by  F.  B.  Jensen  and  W.  A.  Kuperman  [Ref.  9],  with 
predictions  that  satisfactorily  agreed  with  the  experimental 
results  for  small  ray  angles.  With  some  restrictions, 
normal  mode  theory  is  applicable  for  sound  propagation  in 
the  wedge-shaped  fluid  with  a  fast  bottom.  The  parabolic 
eauation  is  good  for  fast  and  slow  bottom,  but  with  the 
restriction  that  horizontal  ray  angles  must  be  less  than 


Another  technique  introduced  to  predict  the  propagation 
of  sound  in  the  wedge  is  the  method  of  images.  This  method 
was  derived  from  the  simplest  case;  a  monofrequency  Doint 
source  in  a  homogeneous  ocean  with  parallel  boundaries.  The 
total  pressure  is  the  sum  of  an  infinite  number  of  spherical 
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the  speed  of  sound  is  constant,  and  that  the  interface  be¬ 
tween  the  fluids  and  the  surface  is  smooth. 

In  both  WEDGE  and  XSLOPE,  all  distances  are  scaled  in 
units  of  the  "dump  distance."  A  dump  distance  X,  as  stated 
in  Reference  10,  is  the  distance  from  the  apex  measured 
alonq  the  interface  at  which  the  lowest  mode  attains  cutoff. 
If  the  wedge  angle  is  8  (Figure  1.1c),  then 


k  ^  sin9c  tan8 


9C  =  arccos ( c^/c2 ) 


(1.1) 


(1.2) 


where  k^  is  the  wave  number  in  the  wedge  and  9C  is  the  crit¬ 
ical  grazing  angle  for  reflection  of  sound  from  the  bottom. 
For  3  <  <  1 


X  =  H  tan  6 


(1.3) 


This  scaling  distance  negates  the  necessity  of  specifying 
f  reauency . 


C.  COMPUTER  PROGRAM  DSLOW 

At  the  start  of  the  work  reported  in  this  thesis,  a  com¬ 
puter  program  was  obtained  [Ref.  13],  which  is  an  extension 
of  the  WEDGE  and  XSLOPE  for  downslope  configuration  with  a 
slow  bottom.  The  computer  model,  DSLOW,  developed  to  run  on 
a  desktop  computer  (Wang  2000),  uses  the  method  of  images  to 
predict  the  pressure  amplitude  and  phase  anywhere 
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within  the  wedge  fluid  overlying  a  slow  bottom  in  a  cross¬ 
slope  configuration.  A  geometrical  picture  of  this  configu¬ 
ration  is  shown  in  Figure  1.1c. 

riathemat ical ly ,  the  model  used  in  WEDGE  and  XSLOPE  is 
applicable  in  any  condition.  But  consideration  must  be 
given  for  nakinq  it  work  for  a  slow  bottom.  In  the  case  of 
a  fast  bottom,  the  dump  distance  has  a  physical  meaning. 
The  dump  distance  is  expressed  as  a  function  of  the  critical 
angle.  The  critical  angle  is  equal  to  arccos  (C]_/C2>.  In 
the  case  of  slow  bottom,  ci/c2  is  greater  than  1,  thus 
arccos  (ci/ct)  is  invalid;  therefore,  so  is  the  dump 
distance.  To  facilitate  the  scaling  factor,  a  "characteris¬ 
tic  distance"  or  "scaling  distance"  is  introduced.  We  need 
the  scaling  distance  because,  with  this  distance,  our  model 
will  be  independent  of  frequency  as  in  the  fast  bottom  case. 
There  is  also  the  hope  that  the  use  of  a  scaling  distance 
will  allow  systematic  observation  of  the  pressure  field. 
This  scale  distance  X0  is  the  distance  measured  along  the 
interface  from  the  apex  to  the  point  where  the  lowest  mode 
would  attain  cutoff  if  the  fluids  in  the  wedge  and  in  the 
bottom  were  to  be  interchanged.  The  characteristic  distance 
is  defined  by  the  following  equation: 


9q  ta 
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where  80  =  arccos ( c2/c ^ )  and  K2  =  w/c2  is  the  wave  number 
in  the  bottom. 

The  following  tens  will  he  used  throughout  (see  Figure 

1.1): 

8  =  wedge  angle 

=  distance  of  the  source  from  the  aoex  in  units  of 
*0 

? =  distance  of  the  receiver  from  the  anex  in  units  of 
^o 

y  =  angle  of  elevation  of  the  source  above  the  bottom 

6  =  ancle  of  elevation  of  the  receiver  above  the 

bottom 

Yq  =  distance  between  the  projection  of  the  source  and 
receiver  on  the  shore  line,  scaled  by  X0 

p^/p2  is  the  ratio  between  the  density  of  the  fluid  in 
the  wedge  ( pj )  and  the  density  of  the  fluid  in  the 

bottom  ( p2  ) 

c ^ /c 2  is  the  ratio  between  the  speed  of  sound  in  the 
wedge  (c^)  and  the  speed  of  sound  in  the  bottom 
( c 2 ) ♦  ^  fast  bottom  occurs  when  c2  >  c^;  a  slow 

bottom  occurs  when  c2  <  ci 

The  mjrnose  of  this  research  is  the  following: 

1.  To  transfer,  test,  and  evaluate  DSLOW  program  cn  the 
IBM  3300; 

2.  To  obtain  numerical  and  graphical  output  for  a 
number  of  cases;  and 

3.  To  attempt  to  develop  plausible  explanations  for  anv 
significant  features  observed. 
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•a.  upslope  configuration 


II.  THEORY 


A.  GENERAL  VIEW  OF  A  WEDGE  PRESSURE  DISTRIBUTION  IN  THE 

DOWNSLOPE  CONFIGURATION 

A  general  picture  of  the  sound  energy  propagation  within 
the  wedge  in  downslope  direction  is  given  in  Figure  2.1.  If 
a  sound  source  is  placed  at  point  S,  ray  1  will  reach  the 
surface  at  point  P  with  an  incident  angle  a  with  respect  to 
the  normal  to  the  surface  at  this  point.  This  ray  is  re¬ 
flected  by  the  surface  at  the  same  angle  but  with  the  phase 
180°  different.  (On  the  surface,  sound  pressure  is  zero 
everywhere.)  The  reflected  ray  reaches  the  bottom  with  an 
incident  angle  8  +  a  .  At  great  enough  distance,  ray  1 
never  reaches  the  bottom  again.  This  ray  does  not  contri¬ 
bute  to  a  sound  pressure  field  at  the  bottom.  The  pressure 
at  the  bottom  should  be  very  small  according  to  the  ray 
theory  argument. 

Using  these  ray-tracing  methods,  an  estimated  profile  of 
the  pressure  amplitude  versus  the  receiver  depth  can  be 
made.  When  the  source  and  the  receiver  are  placed  near  the 
apex,  the  pressure  amplitude  is  zero  at  the  surface,  a  maxi¬ 
mum  somewhere  within  the  wedge,  and  greater  than  zero  at  the 
bottom.  In  the  case  where  the  source  is  at  a  far  distance, 
the  pressure  amplitude  is  equal  to  zero  at  the  surface,  a 
maximum  somewhere  within  the  wedge,  and  zero  at  the  bottom. 


Ray  tracing  will  only  give  a  rough  approximation,  not  an 
exact  solution,  but  ray  tracing  may  be  used  as  a  guide.  The 
method  of  images  calculates  the  exact  pressure  amplitude  at 
each  point  within  the  wedge  subject  only  to  the  assumption 
inherent  in  using  the  plane-wave  Rayleigh  reflection 
coef f icients . 


B.  SOUND  PRESSURE  AT  A  POINT  IN  THE  WEDGE  DOWNSLOPE 

PREDICTED  BY  THE  METHOD  OF  IMAGES 

Let  the  source  be  a  scaled  distance  Rj  from  the  apex  and 
at  an  angle  of  y  measured  from  the  bottom  of  the  wedge.  Let 
the  receiver  be  a  scaled  distance  R2  from  the  apex  and  at  an 
angle  6  measured  from  the  bottom. 

Using  Figure  2.2,  let  the  upper  half  family  of  images  be 
n  =  1,2, 3, 4,...  and  the  lower  half  family  be  n'  = 
1,2,3,4....  Calculating  the  field  resulting  from  source  and 
images  proceeds  along  the  lines  developed  in  [Ref.  14].  If 
$n  is  the  angle  formed  at  the  apex  between  the  nt^  image  of 
the  source  and  the  receiver,  then 
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or  more  generally 

4>n=  (n  +  l)8-6-Y  for  n  odd 
<|>n  =  n0  -  5  +  y  for  n  even 

Which  can  be  reduced  to: 


=  {n+ ( 1/2 ) [l-(-l)n] } 6  +  ( - 1 ) n y  -  6 


(2.1  ) 


or 


$n  =  2  INT [ 


n  +  1 


] 8  +  (-l)nY  -  6 


(2.2) 


where  INT[  ]  denotes  the  largest  integer  which  is  equal  to, 
or  smaller  than  the  argument.  Using  the  same  method  for  the 
member  n'  of  the  lower  family  of  images  we  obtain: 


=  (n+(l/2)[l-(-l)n] }8  +  ( -1 ) nY  +  6 


(2.3) 


or 


♦nt  =  2  INT ( 


n  +  1 


)8  +  (~l)nY  +  5 


(2.4) 


Using  the  geometry  of  Fiqures  2.2  and  2.4,  the  distance 
between  the  nfc^  and  n'^  images  to  the  receiver  is 
respectively 


:n  =  /  +  R2 ^  -  2R^  R2COS  4> n 


2.5: 


and 


rn  1  =  /r"j~2  +  -  2Rj  R2C0S  $n 


(2.6) 
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respectively  are 


9no  =  arc  tan  [ 


sm  *n 


R2/p1  “  ccs  4»n 


(  2.7  ) 


an.' 


.  ?  i  n  4>n  i 

9ni„  =  arc  tan  - — - LI - 

n  °  *2/ll  ~  cos  ♦n' 


(2.8) 


Define  9nrT,  and  9pin  as  the  anales  of  incidence  for  the 
bounces  from  the  bottom  for  the  n  and  n'  image 
respectively;  m  =  1,2,3....  (The  0th  bounce  is  the  last 
one  before  reaching  the  receiver.)  The  geometry  of  Figures 
2.2  and  2.3  nive  9nP  as  follows: 

e21  =  920  '  26  “  6 

03 2  =  930  ~  43  ”  6 

041  =  © 4 ij  ~  2  8  ~  6 

95 2  =  95U  “  40  "  5 


The  general  expression  is 


9nm  9no  ~  2rn9  “  ^ 

Using  the  same  method 

9n  '  m  9n  '  o  2m  9  +  $ 

The  maximum  number  of  bottom  bounces  of  the  nth  and  n,th 


image  is 


IST  Un/2B]  = 


max 
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INT  [  4»n  »  /2  8] 
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The  maximum  number  of  images  is 


nm*v  =  N  =  INT  [tt/8! 


The  reflection  coefficients  for  the  nfc^  and  n  ’  i 
for  a  plane  wave  are: 


mages 


where 
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T2^2 


P2c2 


ROn'm)  = 


J2^2_ 

9 1.9.1 

p2c2 


-  Yn'i 


+  'I'n'i 


( 2.10  ) 


/  1  -  (C]/c?)2COS2  9~ 

sin  9nm 


(2.11  ) 


1  -  (C]/Cy)^CQ5Z 
sin  9n ' m 


(2.12) 


The  contribution  from  the  upper  family  of  images  is 


Pu  =  l  -L-  exp(-jkrn){-l)INT[n+1)/2l  nRnm  (2.13) 
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DS LOW  PROGRAM  IMPLEMENTATION 


A.  PROGRAMS  FEATURES 

Since  the  mainframe  graphics  computer  was  available,  the 
DISSPLA  graphical  program  was  used.  The  only  programming 
language  compatible  with  DISSPLA  is  FORTRAN.  The  numerical 
and  graphical  output  is  provided  by  this  program.  To  give 
the  pressure  amplitude  versus  received  angle  graphs,  two- 
dimensional  plotting  is  used. 

The  program  DSLOW  is  run  by  placing  the  point  source 
anywhere  in  the  wedge  and  then  placing  the  receiver  at  a 
distance  downslope  from  the  source.  The  receiver  position 
was  varied  from  zero  degrees  at  the  bottom  to  6  at  the 
surface.  High  resolution  plotting  was  achieved  by  dividing 
the  y-axis  (received  angle)  into  two  regions.  The  first 
region  covers  the  receiver  angles  from  zero  to  1/5  of  the 
wedge  angle.  In  this  region  A5  is  equal  to  6/100.  The 
second  region  covers  the  remaining  wedge  angle  with  A6  equal 
to  6/10.  This  method  provides  29  predictions  of  the  pres¬ 
sure  amplitude.  Another  method  of  plotting  carried  out  was 
in  the  region  of  6  >  6/2,  A5  =  8/10,  and  in  the  region  of  <5 
<  6/2,  A<5  =  6/100.  This  method  provides  54  points  to  be 


B.  NORMALIZATION 

The  main  goal  of  this  research  was  to  investigate  the 
profile  of  the  pressure  amplitude  as  a  function  of  a  number 
of  variables.  An  example  of  the  numerical  values  of  the 
pressure  amplitude,  the  normalized  pressure  amplitude,  and 
the  phase  at  each  receiver  position  is  displayed  in  Appendix 

C.  The  sound  pressure  becomes  smaller  as  the  receiver  is 
moved  away  from  the  source.  If  the  pressure  amplitude  were 
clotted  directly,  it  would  be  difficult  to  compare  the 
curves  at  near  distances  to  the  curves  at  far  distances 
since  at  the  near  distances  the  pressure  amplitude  is  much 
greater  than  the  pressure  amplitude  in  far  distance.  Thus, 
a  normalized  oressure  amplitude  is  needed.  The  normalized 
pressure  is  obtained  as  follows:  (see  Figure  3.1) 

We  know  that  the  sound  pressure  at  the  surface  is  zero 
and  that  the  sound  pressure  is  a  small  number  greater  than 
zero  at  a  point  near  the  surface.  The  first  non-zero  value 
of  pressure  P]_  is  at  the  receiver  angle,  6^  =  9(5/10.  We  use 
this  first  calculated  non-zero  pressure  amplitude  as  the 
normalization  unit.  The  normalized  pressure  is 

PN  =  P( 6)/( Pl )  (3.1) 
where  P(5)  is  the  pressure  at  any  point  within  the  wedge. 
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C.  PROCEDURE 

Figures  3.2  through  3.7  represent  the  results  when  the 
receiver  distance  and  source  angle  are  fixed  and  the  source 
distance  and  receiver  angle  are  varied.  Figures  3.8  through 
3.13  represent  the  results  when  the  source  distance  and 
angle  are  fixed  and  the  received  distance  and  angle  are 
varied.  These  cases  will  be  the  foundation  of  our 
subsequent  discussions. 

The  solid  lines  indicate  the  fitted  curve  and  the  dots 
indicate  some  values  of  the  normalized  pressure  amplitude. 
In  DS LOW ,  the  dot  appears  at  each  third  datum. 


D.  PROGRAM  IMPROVEMENT 

DS  LOW  was  desiqned  to  provide  three-dimensional  graphs. 
For  example,  the  x-axis  represents  the  scaled  source 
distance,  the  y-axis  represents  the  scaled  received 
distance,  and  z-axis  represents  the  normalized  pressure 
amplitude.  To  simplify  the  presentation,  only  two- 
dimensional  graphs  were  presented  with  the  x-axis  the  nor¬ 
malized  nressure  amplitude  and  the  y-axis  the  receiver  ancle 
5.  a 1 1  curves  are  presented  with  the  data  fitted  with  a 
cubic  spline. 

The  DSLOW  program  was  executed  to  obtain  numerical 
results  of  the  phase  angle,  the  pressure  amplitude,  and  the 
normalized  pressure  amplitude  at  each  receiver  position. 
The  first  run  used  double  Precision  for  accuracy. 
Difficulties  were  encountered  when  the  DISSPLA  subprogram 


3 


was  attached  for  making  the  graphical  output.  When  double 
precision  and  DISSPLA  were  not  successful,  the  single  pre¬ 
cision  was  used,  resulting  in  round-off  error.  (See  Figure 
3.8  at  R2  =  10.0.) 
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Figure  3.2  Graphs  of  receiver  angle  6  versus  pressure 
amplitude  with  F*2  fixed,  varied 
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Figure  3.3  Graphs  of  receiver  angle  6  versus  press 
amplitude  with  f*2  varied 


NORMALIZED  PRESSURE  AMPLITUDE 


RECEIVER  ANGLE  VS.  PRESSURE 


^  § 

H  I"  B 

11^1 

»  *c8 


o  o 
w  s 


b  6  _  8 

5  5  8° 

BfceU? 

as^ls 

c*3  u  ^  ** 

^s§ 

BS^Ss 

So  55  o:  c/i  ^ 


o 

■*a 


§2 

~  D 


o  ft- 
•cd  ft. 


°§ 


0'9  0>  02 

(03d)3TDNV  H3AI3D3H 


Figure  3.5  Graphs  of  receiver  angle  6  versus  pressure 
amplitude  with  R2  fixed,  R}  varied 
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Figure  3.6  Graphs  of  receiver  angle  6  versus  pressure 
amplitude  with  R2  fixed,  R^  varied 
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Figure  3.7  Graphs  of  receiver  angle  6  versus  pressure 
amplitude  with  R2  fixed,  Rj  varied 
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Figure  3.8  Graphs  of  receiver  angle  5  versus  pressure 
amplitude  with  R2  fixed,  Rj  varied 
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Figure  3.10  Graphs  of  receiver  angle  6  versus  pressure 
amplitude  with  R2  fixed,  varied 
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IV. 


DISCUSSIONS 


A.  GRAPHICAL  OUTPUT 

The  qranhs  of  normalized  pressure  amplitude  as  a  func¬ 
tion  of  receiver  anale  5  were  investigated  for  various 
source  distances  and  receiver  distances  R2  (Fiq.  3.2  to 
3.13)  while  the  other  parameters  are  held  contstant.  For  a 
given  wedge  angle  3  and  sufficiently  small  source  distance 
(Fiq.  3.8  and  3.12),  at  all  receiver  distances,  the  pressure 
increases  uniformly  towards  the  bottom.  For  greater  source 
distances,  (Fiq.  3.9  -3.11)  the  pressure  attains  a  maximum 
within  the  wedqe  for  all  receiver  distance. 

As  the  receiver  distance  is  increased  (Fiq.  3.9),  a 
pressure  minimum  develops  between  the  maximum  and  the 
bottom.  An  important  property  of  the  curves  of  pressure 
versus  receiver  angle  when  there  is  a  maximum  and  minimum  is 
that,  at  a  specific  receiver  distance,  the  pressure  above 
the  mjriipum  can  be  extrapolated  to  zero  pressure  on  the 
bottom.  (See  Fig.  3.10  with  R2  =  32.)  This  receiver  dis¬ 
tance  is  called  the  "transition  point."  So  far,  we  do  not 
know  the  properties  of  the  transition  point.  We  use  the 
transition  point  for  indicating  the  behavior  of  the  curves 
when  the  parameter  involved  is  varied.  The  transition  point 
appeared  twice  in  some  cases,  but  in  the  following 
discussions  the  first  transition  point  is  the  only  point  we 
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will  be  concerned  with. 


(See  Fig.  4.1  for  transitions  cor¬ 


respond  to  R2  =  4.6  and  6.4.) 

B.  GRAPHS  CLASSIFICATION 

The  development  of  curves  with  the  source  distance  (R^) 
and  the  receiver  distance  (R2)  as  variables  was  observed. 
As  Ri  or  R2  are  varied  the  curve  changes  from  a  linear  curve 
to  a  curve  with  an  observable  minimum  (Fig.  3.9,  R2  =  5.0) 
and  finally  to  a  curve  without  a  minimum  (Fig.  3.9,  R2  = 
9.0).  Three  different  types  of  curves  resulted  from  the 
series  of  two-dimensional  plotting.  They  are  described 
below: 

1.  Type  1  Curves 

Type  1  curves  (Fig.  4.2)  are  those  where  the  sound 
pressure  is  equal  to  zero  at  the  surface  and  maximum  at  the 
bottom  and  is  almost  linearly  dependent  on  depth.  These 
curves  are  most  pronounced  when  the  source  distance  is  much 
smaller  than  the  characteristic  distance.  The  closer  the 
source  is  to  the  characteristic  distance,  the  more  nonlinear 
the  curves  (Figs.  3.3  and  3.4). 

2 .  Type  2  Curves 

Type  2  curves  (Fig.  4.3)  are  those  where  the  sound 
pressure  is  zero  at  the  surface,  maximum  somewhere  between 
the  surface  and  the  bottom  with  no  minimum.  These  types  of 
curves  are  generated  when  the  source  is  placed  at  a  point 
much  greater  than  the  characteristic  distance.  Type  2 
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curves  indicate  that  the  sound  energy  in  the  wedge  is  well 
colligated  and  that  reflection  is  negligable. 

3.  Type  3  Curves 

Type  3  curves  (Fig.  4.4)  are  those  that  have  a  min¬ 
imum  pressure.  These  curves  occur  when  the  source  is  a  dis¬ 
tance  slightly  greater  than,  or  less  than,  the  charac¬ 
teristic  distance.  Tables  1,  2,  and  3  of  Appendix  D  show 
the  receiver  positions  at  the  first  transition  points. 
Three  different  values  of  8,  two  different  values  of  p[/P2> 
and  two  different  values  of  c^/c 2  were  used  in  making  these 
tables.  The  transition  point  did  not  occur  when  8  =  15°, 
Pl/?2  ~  0.90,  c^/c2  -  1.10.  An  explanation  can  be  offered 

using  the  fact  that  for  these  particular  sound-speed  and 
density  ratios  an  angle  of  intromission  exists  [Ref.  15]. 
Since  the  angle  of  intromission  is  the  grazing  angle  at 
which  the  sound  energy  is  completely  transmitted  into  the 
slow  bottom,  it  is  plausible  that  no  transition  point 
occurs. 
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C.  TRANSITION  POINT 

By  varying  the  wedge  anale  8  in  small  increments  AS  - 
u.5°  starting  with  8  =  5°,  and  ending  at  8  =  7°,  it  was 

found  that  transition  occurs  for  source  distances  within  the 
range  from  1.0  to  1.5. 

por  <  l.o,  no  transition  point  was  observed;  the 
curves  are  the  Type  1.  For  1.0  <  Ri  <  1.5,  the  evolution  of 


Figure  4.3  Type  2  curves,  indicating  a  well 
sound  field  as  the  source  away  f 
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D.  PARAMETER  VARIATIONS 

Variation  o£  parameters  was  done  by  changing  one 
oaraneter  of  interest  while  all  others  were  held  constant, 
for  fixed  source  and  receiver  distances,  and  plotting  the 
receiver  angle  versus  normalized  pressure  amplitude. 

The  parameters  3,  oi'o2>  and  c^-^,  were  hell  constant 
and  the  pressure  amplitude  was  clotted  for  various  ?,  Rj/ 
and  R2 ♦  The  y’s  are  set  at  3/4,  3/2,  and  36/4.  Variations 
in  the  shore  distance  (Yc)  can  be  made  because  the  program 
is  available,  but  to  simplify  the  investigation,  Y0  was  set 
equal  to  zero  for  all  plots  (Fig  4.2  is  included  as  an 
example  for  Y0  *  0). 

1.  Variations  of  6 

Initially,  the  values  of  3  investigated  were:  6°, 
10°,  and  15°.  The  major  effect  created  by  altering  the 
value  of  8  is  that,  for  the  sane  values  of  Rj  ,  p^./p2»  and 
C\:'C2'  t-^ie  seller  8,  the  shorter  the  transition  point  (see 
Figs.  3.9  and  3.10,  Tables  1,  2,  and  3  of  Appendix  0). 

2 .  Var iat ions  of  y 

The  variations  of  y  from  y  =  6/4  (the  source  is 
placed  near  the  bottom)  to  y  =  38/4  (the  source  is  placed 
near  the  surface)  are  presented  in  Tables  1,  2,  and  3  indi¬ 
cated  that  the  greater  y  the  shorter  the  transition  point. 
It  is  not  always  true,  for  instance  in  Table  1  at  8  =  6°,  R^ 
=  1.50,  the  greater  y  the  longer  transition  point,  for  the 


rises.  (See  Appendix  D. ) 

3.  Variations  of  c\/c 2  and  p^/p2 

Variations  of  the  acoustical  parameters  c\/c2  anc^ 
Pl/P2  were  done,  but  did  not  give  a  significant  variation  of 
the  sound  pressure  profile. 
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V.  CONCLUSIONS  ViD  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  progran  DSLOW  gives  good  plots  representative  of  the 
sound  energy  distribution  within  the  wedge.  The  sound  energy 
can  be  well  collimated  by  the  wedge.  This  phenomenon  is 
strongly  affected  by  the  source  position.  At  a  source 
position  close  enouoh  to  the  apex,  sound  energy  is  distri¬ 
buted  linearly  with  respect  to  the  depth.  As  the  source 
moved  away  from  the  apex,  the  distribution  of  sound  energy 
becomes  more  complex.  Sometimes  a  minimum  is  found;  this 
minimum  may  be  caused  by  the  presence  of  sound  energy 
reflected  by  the  bottom. 

The  source  position  plays  a  maior  role  in  forming  the 
pressure  distribution  profile.  The  pressure  distribution  is 
also  very  sensitive  to  the  parameter  variation  at  small 
source  distances,  but  it  becomes  insensitive  at  large  source 
distance.  The  characteristic  distance  must  have  physical 
meanings  rather  than  just  an  arbitrary  number,  because  when 
the  source  distance  in  proximity  to  the  characteristic 
listance,  the  model  is  most  sensitive. 

The  model  is  restricted  when  the  single  precision  mode 


generates  round-off  error  and  rough  curves  which  do  not 
allow  for  accurate  analysis. 


RECOMMENDATIONS 


1.  Sinale  precision  Produces  aocd  results,  but  failed 
in  some  cases.  Double  precision  would  improve  the 
procram,  but  increase  the  execution  time.  This  must 
be  done  by  runninq  the  nrocra^  in  double  precision, 
and  accumulatinq  the  result  in  sinale  precision 
before  olottina  the  data  by  DISSPLA. 

2.  It  is  sunoested  that  the  program  be  run  using  more 
realistic  parameters  and  observing  the  effects  on 
the  characteristic  distance  and  transition  point. 

3.  Further  studv  validating  DSLOW  in  comparison  with 
experimental  results  is  suggested. 


APPENDIX  A 

DSLOW  ALGORITHM 

The  pressure  amplitude  calculation 

N,  =  I NT[  1  SO  P) 

(eqn  A.l ) 

K 

K  i  X  = - - - 

2tanP  tan[arccos(c->  Cj)] 

(eqn  A. 2) 

AL  =  a  K->  =  0.0001  (constant) 

(eqn  A. 3) 

d2">o^r12-R22 

(eqn  A. 4) 

R.-2R|R, 

(eqn  A. 5) 

Rg  =  v  Dt  -  cos[(N-l)P-t-y  -  6| 

(eqn  A. 6) 

Rg  =  v  Dt  -  R3  cos((N-l  )P  +  y  6] 

(eqn  A. 7) 

S2  =  (.i)INT(Ni2) 

(eqn  A. 8) 

\Vj-2AKcj  c2)2 

(eqn  A. 9) 

SI  =  1  { R j  sin[(N-l)P  +  y)  -  2(INT((N-1)  2)P 

(eqn  A.  10) 

+  R2  sin(2INT((N-l).2)P-6]J!  Rg 

Cl  =  v  ( 1-SI") 

(eqn  A.  1 1 ) 

* 

T  =  SI  Dj 

(eqn  A. 12) 

Wo  =  (‘c2't'cl  c2) 

(eqn  A.  13) 

Y  =  VW0‘  +  \vy 


(eqn  A. 14) 


Z=IW0I 

(eqn  A.  15) 

V^v'lY  +  ^l  2 

(eqn  A. 16) 

Y2  =  V(V-VY0),2 

(eqn  A. 17) 

.  _  Tj-Y, 

1  (T-Y2r  +  Y^ 

(eqn  A. IS) 

v, 

- i - r~ 

-  (T-Y2r-Yj- 

(eqn  A. 19) 

,  _  (T2-V2)2- Y,2 

3  (T-Y2r-Yr 

(eqn  A. 20) 

2V,T 

ZA  = - , - r5 - 1 

6  (TMSr-Y^ 

(eqn  A. 21) 

N 

1 

Pl=  I(-l)INTfN  2)(Z5  cos(R8nK1X)  +  Z6  sm(R8nK,X}  Rgn 
n  =  1 

(eqn  A. 22) 

N 

1 

P2  =  I(-D1NT(N  2){-Z5  sin(R8nK1X)  +  Z6cos(R8nK,X}.R8n 

n  =  1 

(eqn  A. 23) 

N 

1 

P3=  V  (.i)IM(N  2){Zj  C0S(RgnK|X)  +  Z£  Sm(R^nKjX}  R9n 
n  =  1 

(eqn  A. 24) 

m 


» 


feqn  A. 25) 


V  (.1)1NT(N:2){.Z5  sin(R9nK,X)  +  Z6  cos(R9nK,X)  R9n 
=  1 


Pi  +  P2 


P6  =  P3  +  P4 


Ptot=Rlv/p?7P7 


(eqn  A. 26) 


(eqn  A. 27) 


(eqn  A. 28) 


APPENDIX  B 

DSLOW  PROGRAM 


************************************************************************ 
**  THIS  PROGRAM  CALLED  DSLOW  IS  CALCULATING  THE  SOUND  PRESSURE  ** 

**  WITHIN  THE  WEDGE  OVERLYING  SLOW  BOTTOM  FLUID  AND  DOWN-SLOPE  ** 

*****************************$kS^**k$*********************************** 


INTEGEF 
,  REAL"? 


REAL *4  TOQ 


A,I,I1,M.N,S1,S2.N1, 


1,TQS2,IM3 


******«K*K*KtK*K^RK*W******^***************'**llr^*****yr******  ****★***★★★* 

******^»T«yX*KK^^XKK^^^********** ******* ********************************** 
C  B  =  WEDGE  ANGLE  (DEG) v 
C  G  =  SOURCE  ANGLE  7dEG) 

C  D  =  RECEIVER  ANGLE  7DEG) 

C  Nl=  fr  OF  IMAGE  POINTS 

C  R1  =  SOURCE  DISTANCE  (IN  CHARACTERISTIC  DISTANCES) 

C  R2  =  RECEIVER  DISTANCE  (lN  CHARACTERISTIC  DISTANCES) 

C  YO  =  APEX  DISTANCE  (IN  CHARACTERISTIC  DISTANCES) 

C  D1  =  RHO  1/RHO  2 

c*****t^i . ... 


C  AL  =  ALPRA/K2 

y*****^****^*y y*xx yy jyy x*^y y k*w5 ************************************** 
****** 5 Ky**y?KK*§**yjfx^i§*y ***************************************** 

S  IS  :  i:88 

5  gVii.r1-6 

B  =  10.0 

S  8 :  m 

--  i ;  J y% 

C  R1  =  6.95 

R1  =  1.1Q 


:  m 


K4  =  3.U 

yo  =  o.g 

D1  =  0.8 

li  l  4:4 


■k-k -kick -knit  ★★★***★*★  *★★★★★★★★★★★★*★  ***★★*★★ 

★  *  *  ★  ★  -k  k  «  k  -k 

C2  =  CC**2 

3  D2  =  (Y2*Y2^+(RI*R1)+(R2*R2) 

c  iL*  Pi/(2*TAN(ACOS(l/CC))*TAN(B)) 


=  PI/ (2*TAN(AC 
2  =  TMHBl  ,  , 
n  =  Acosfi/cc 
52  =  tanTtooiT 

=  PI^TQQj  5 


60 


.FORMAT 

WRITE11 

« \mT, 

FORMAT 
WRITE  i 
^FORMAT i 

WRITE  i 
FORMAT  I 
WRITE  I 
FORMAT ( 
WRITE 
FORMAT 
WRITE 
FORMAT] 
*  1  PHi 

WRITE  < 
FORMAT 
WRITE 
DO  10 


J|E"|r  R9  .J-.'THETMH)  ',7X,'IMGE  SR  R8  1  ,7X, 

6 , 24 1 1 

6,2f8lEE$r£LG*T6  ' .F5.2.1X,'  SOURCE  ANGLE  =  '  ,F5 .2) 

'  IhorFdISTANCeI'!  'RECEIVER  DISTANCE*  ',F6.2,1X, 

6  251;  R1  R2  YO 

6§70)/Dl<iCC,AL5‘2'5X' '  C1/C2=' .FS.2.5X. 1  ALPHA/K2*' ,  F8.4) 

:'kix  _=,c?:f12) 


1 , F5 . 2 , IX , 1  SOURCE  ANGLE  =  -,F5.2) 

;  ,F4. 2  IX, 'RECEIVER  DISTANCE*  ',F6.2,1X, 


6,272) 
1*, 'REC, 
SE  ANGLE 

,6,800) 


'  •2X-  'PRES.MIPLIIUDE 


:TE  6,801 
10  M>  0  ,  A 
D  =  M*B/A 
DD=B/A 
V  =2*8/10 
IF  D.LT.V)  C 
IF(D.GE.V)  C 


IF(D.G 
DX  =  D 
DO  15  J 
D  =  D 
DD  =  B 
SI  =  1 
DO  20 

mi 

R^(N 
R9  N 


:  f?jskl!*B/(10‘A) 

111)  =  s^rtJdI-r^coKtiInI+d! 


R|(N)  =  SQRT(D2-R3*C0S(T1{N)-D)) 
FORMA?fl&2  %  p uh?wM+n\) 
CCNTINuI '  3l° }  ** ' 11  (N 5 '  ^  ^  W  ^ 


n  =  o. 

52  f  0. 


’w?°:  !4M™T(N/2)) 

$UP***® 

wo  =  |-C2+7C(I)*C(I))) 

IF(YTLE,Zl  y  =  z 
Y1  =  Ql*s6RT(Y+W0L 

Z1  =  f?^2^2R^y'W^ 

Z2  =  -Yl 

Z12*  CMPLX|51R  Z2J) 

24  =-Z2^Zl*Zl+H*Z2) 

If3:  §®<IP44I) 

Zl|=  CMP|X(Z5RZZ6I) 

kU‘  ,§9?x  ,^5X,  'I  '  , 5X ,  1  ARSIN(S ( I )*T6) '  ,5X,  'E(I) 


:W^)  +  (wl^l 


:F0RMi^ 

WRITE  (6,400)  ^ 


( I )+F( I ) *F( I ) ) 1 , 5X, 1 T6*ATAN(F(I)/E(I) ) 1 ) 


C420  FORMAT 
C  ^WRITE 

-  40  CONTIN 


3X , F6 . 4 , 10X , F7 . 4 ) 
(F(I)/E(D) 


IF(N .LE .2.00)  GOTO  50 

°  2^=  Ip1'  11 

It  =  il*Z3-22*Z4 
CONTINUE  I  21  W*z! 

z\  =  p 

h  =  ?4os?fSN) 

it  E  iiw-ktt 

PI  =  pj+sf*Z5/R8(N) 
=  P2+S2*Zo/R8(N 5 

mj uw  i*p 


FORMAT ( 
WRITE  ( 
FORMAT ( 
WRITE  ( 
^ FORMAT ( 

WRITE  ( 

DO  60  I 


:E(REFL)=  ' ,3X, 'IM(REFL)* 


*slh%  IS*!i?fo£!i4> 

,^'ARCSI^S(lt)*T6=*,3X.' 


5X  ';aKI2 
IS(li '=:ABS 

T^SClf?!) 

WO  =  -C2+C 

z :  Blftt 

IF  Y.LE  zY 

B :  84*B 


E (I ) 1 ,5X, ’F(I) 1 , 5X, 'EF1 


w^wS)+|wi*wi)) 

Y  =  Z 
RT (Y+WO) 
QRT(Y-WO) 


T+Y2 

Yl^ 

Z1*Z3-Z 

Z1*Z4+Z 


AN1-  ASIN(S ( I 


LO  FORMAT  I 
WRITE  l 
60  CONTINl 


nm  —  noiiNijiiJ  y  i  q 

E2: M  §  5  /!  i]Wf>*F<I» 

1 S ; 4  ¥ ;  bi  f  e?5  j  r  I  { ! ;  Sgf  f  In!  - 5X  • F  • 4  ■ sx  - 


F7.4) 


62 


C  700 

c 

C710 

C 

-  30 


'11° 

-10 


2  =  26^ 

:  ”^?r 

=  21^Z3-l2*Z4 

pi  =  pi+i^p/Ri(N) 

" 2  =  P2+S2*Z6/R9 {N i 

PATH  NO=  * ,  3X , ' RE ( REFL )  = 


VtSh  HIM6-* 


FORMAT 
WRITE 
FORMAT 
WRITE 
CCNTINL 

^7 ( R^=n*T6 

PZ  K  =SQRT(P1*P1+P2*P2)*R1 
AT (K2=Af ANiPZ/Plj 
WRITE  (6^310;  K,D 
FORMAT 
CONTINUE 
CONTINUE 

r.n  :i  T  —  1  £9 


, 3X , ' IM(REFL)=  '  ) 


Xp=FZ(£ll ' 
pn(lt=??(;)/kp. 

VRI 


rfi1 

C 

-111 


knk  n-k 
kx 

kkkkk 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


WRITE 
FORMAT 
CONTINUE 
Y0=YQ+0.5 
R2=R2*0.50 
CONTINUE 
STOP 


A  SUBPROGRAM .FOR .PLOTTING . BY . TEKo 1 80R . SWERPA 


** 

******************** 


CALL  MEDBUF 
CALL  TEK618 
CALL  COHPRS 
CALL  NOBRDR  ,  „  . 
CALL  PAGeTiS .  ,  12  . ) 
CALL  AREA2D(ll, ,9.) 
CALL  HEIGHT [ . 2] 

CALL  . .  ‘ 

CALT 

CALL  ’/TICKS  ( 

CALL  XTICKSL.  „ 
CALL  GRAF ( 0 . ,0.10 
CALL  GRAF (0. ,10 
CALL  GRAF (0. ,3 
CALL  DOT 
CALL  GRID (2 , 2 ) 

CALL  HEADIN  1  RE C. 
CALL  ME S SAG 
CALL  MESSAG 
CALL  MESSAG 
CALL  MESSAG 
CALL  MESSAG 
CALL  MESSAG 


-  Jjp  j  /-rj-p  r  2  ) 

L  XNAME( 'NORMALIZED  PRESSURE  AMPLITUDES ', 29 ) 

L  YNAMEl 'RECEIVER  ANGLE(DEG) $  1 , 19) 

^EC . ANGLE  VS.  PRESSURES ', -100 , 1 . 8 , 1 ) 
WEDGE  ANGLE:  S'  100.3. .7.) 

15. OS'  .100, 'ABUt'  , 'ABUT' V 


6.5 

BUT 


> 


CALL  MESSAG 
CALL  MESSAG 
CALL  MESSAG 
CALL  MESSAG 
CALL  MESSAG 
CALL  MESSAG 
CALL  MESSAG 
CALL  MESSAG 
CALL  ME 
CALL  ME 
CALL  ME 
CALL  ME 
CALL  ME 
CALL  ME 

CALL  ME _  , 

CALL  MESSAG 
CALL  MESSAG 
CALL  MESSAG 
CALL  MESSAG 


1/R ft02:  S'  100.8. 

0.80  S'  ioo. 'Abut' , 

C17C2:  S'  .l0O,8.  ,6.0; 

I. 50S'  ,100 'ABUT1 'ABUT' ) 
SOURCE  ANGLE:  f '100.8 ..5,5) 

II. 25$'  ,100, ' ABUt 1 , 'ABUfr 


1  SOURC 
0.25 


ABUT' , 'ABUT' 
'ABUT' , 'ABUT' 
'ABUT'  'ABUT' 
:  S'  100.8..5 
'ABUt' , 'ABUt' 


■SSAGi 
I5SAG I 
SSAGI 
SSAGI 
•SSAG 

•ssagI 

■SSAGi 


1.10  S' ; 100. 'ABUT' . 'ABUT' 
1.6Q$r, IOO, 'ABUT' , 'ABUT1 ) 
0.95  S'  IOO. 'ABUT' . 'ABUT' 
'REC.DISt .INlT:  S' .100  ' 

'0.56  S'  ,100,  'ABUT' , 'A _ 

'  1 . 65  S'  100,  'ABUT'  ,  'ABUT' 
'3.00  S'  100,  'ABUT'  ,  'ABUT' 
'3.00  S '  100 , ' ABUT '  ,  1  ABUT ' 
'REC.DISt. INtRMT  :  S', 100. 

'IW5SP '  i^00, ' ABUT’  ■  'Abut 


Si, 4-5) 


'.ABUT , 

■abut' 


;  |HggD|$T.  jglT>3.li;  100.  s 


'ABUT1 
'0.8.  . 
ABUT' 


.  ,4.0) 
•  5) 


63 


.v.y 


■VIWriurf‘V«  »«< 


APPENDIX  C 

NUMERICAL  RESULTS  OF  DSLOW 


WEDGE  ANGLE  =  10. 0C 

SOURCE  d!STANC|=0.75 
RH31/RH02=  0.80 
K1X  =  19.44 


SOURCE  ANGLE  =  2.50 

RECEIVER  DISTANCE=  9, 
Cl/C2=  1.10 


00  SHORE  DISTANCE=  0, 
ALFHA/K2=  0.0001 


REC.POS  REC. ANGLE  PRES .AMPLITUDE 

i  8:?8  8:88?5 


8:18 


Q .  007175 
0.006965 
0.006763 
0.006563 
0.006374 
0.006 
0.006 
o.oor 
o.oo; 
0.00? 

::§e; 

O.GCL.w. 
0 .  C04§83 
0.004866 
0.004759 
0 .004658 
0.004566 
0.004431 
0.00400- 
0.00391 

o.oo:* 

g  .  OGL  _ . 
0.002981 
0.002149 
0.001121 
0.000000 


PHASE  ANGLE 


-0.732669 

-0.||||4| 

-0 ! 583461 
-0. §45011 
-0.505302 
-0.113330 
0.201979 
0 .41323 
0.5473? 
0.6320] 
0.68086, 
0.69770? 
-0.048596 


NORM. PRESS 


.  _  '227? 
,.497443 

i'.illlu 

!:K98« 

1.000000 

0.000029 


00 


65 


WEDGE  ANGLE  = 
SOURCE  DISTANC 
RH01/RH02*  0 . £ 
XIX  =  19.44 


s-4?*8‘ 


SOURCE  ANGLE  =  2.50 

RECEIVER  DISTANCE*  4.50 
C1/C2*  1.10 


REC.POS  REC. ANGLE  PRES .AMPLITUDE 


6.00 

hi 

10.00 


mm 


u 

0 !  03084€ 
0 .030532 

8:8  S  i 


0.01426S 


0.005003 

0.000000 


SHORE  DISTANCE*  0.00 
ALPHA/K2=  0.0001 


PHASE  ANGLE 


-0.872103 

■0 ! 7^4721 
-0.682756 
•0.660769 
0.273944 


NORM. PRESS 


ii 


62953 
.491957 


6 . 03  So 

I  *  fit? 


1 1 000000 

0.000001 


WEDGE  ANGLE  =  10.00  SOURCE  ANGLE  =  2.50 

?f?clSE?.?iSTANCE=  2'2S  M/R 

REC.POS  REC. ANGLE  PRES .AMPLITUDE 


M/Ril^oo?'00 


9.00 

10.00 


0.220096 

0.219349 

8:11  m 

0.216893 

8:|iiB?8 

0.214095 

0.213090 

0.212046 

oiloilll 

0.208671 

0.194767 

0.176765 

0.154712 

0.128893 

0.099744 

0.068024 

0.034481 

0.000000 


PHASE  ANGLE 

1.084841 

1.092171 

1.099430 

1.106651 

1.113891 

i : iilfli 

1 .1424^ 
1 1 156374 

1:1! 8?l! 

1.176949 

i:il3fl| 

1.197003 

1.203545 

1.210009 

1.216443 

1.222806 

1.231938 

1:111! Ik 

l:tl?l!f 


l:t!l|!| 

■1.243458 


NORM. PRESS 

mm 

6.476104 

6.459502 

6.441518 

!:1  1! 

f:  I  f 

6.237287 

!:?!§!§ 

6.149565 

&:o!lLo 

6.051688 

5.648450 

5.126376 

4.486338 

3-738033 


2.892702 

1.972762 

1.000000 

0.000001 


■  ^  \  V  7.  VV  V 


V.  /  t* 


WEDGE  ANGLE  =  10A00  SOURCE  ANGL^  =  2  1 

I88??Ih8  1  Wj  • ' 

VM  Y  =  19,44  '  ± 


1.13  SHORE  ,DISTANCE=  0.00 
ALPHA/K2=  0.0001 


REC.POS  REC. ANGLE  PRES .AMPLITUDE  PHASE  ANGLE 


5.00 

6.00 

7.00 

S.OG 

9.00 

10.00 


1.900229 

1.923264 

1:111111 
2 1001945 
2.018349 

1:81  m 

2.06002S 

]:pj 

hum 

2.1195S0 

2.122012 

i : iiliio 

2.074631 

1.939167 

1.730041 

1.459312 

mm 

0.‘^00§01 


-0.125102 
-8 1 114742 

:§:i  mi 

-0 .098897 

I# 

-0.061446 

-0.‘ 04^2^1 

-0.034997 

-0.029773 

-0.024576 

-0.019431 

8:8?I|8| 

Q. 112477 
0.144356 


NORM. PRESS 
4.791066 

I  mil 

Mill 

nm 


molt 

344040 

350248 


\\lml 

him 

1.000000 

0.000002 


WEDGE  ANGLE  f  10.00  SOURCE  ANGLE  =  2.50 

SOURCE  DISTANCE^.-) 5  RECEIVER  DISTANCE=  0.56 

RH01/RH02=  0 . 80  Cl/C2=  1.10 

K1X  =  19.44 


SHORE  DISTANCE2  0.< 
ALPHA/K2=  0.0001 


REC.POS  REC. ANGLE  PRES .AMPLITUDE 

1  0.00  3.992624 

i  8:i8 

!  8:18  f:85?§ 

6  0.50  4.093628 

7  0.60  4.10743? 

8  8:18  i:H88? 

10  0.90  4.139806 

11  1.00  4.147007 

12  1.10  4.152471 

13  1.20  4.156213 

1:28  tmn 

if  1:18  t:li!lW 

18  1.70  4.149846 

19  1.80  4.143713 

20  1.90  4.136028 

21  2.00  4.126792 

22  3.00  3.953794 

23  4.00  3.647243 

If  1:88 

h  8:88  i-MVii 

28  9.00  0.728641 

29  10.00  0.00G000 


PHASE  ANGLE 


0.36924? 

0.373416 

0.377565 

0.38i6|| 

0  I  3 §9 8 /b 
0.428945 
0.464082 
0.494513 
0.519774 
0.539597 
0.553331 

•8:il«li 


NORM. PRESS 
5.479548 


1.661543 

,.691426 


1.70684^ 
.707335 
.705572 
.701555 
.695323 
.6^6906 

illllli 

4.425613 
3.71000? 
2.883053 
1.970616 
1 .000000 
0.000000 


-V- '  ".V.  V.V.Wa'.V. 


/ 

✓ 

* 


.’  V  / 


R2  AT  TRANSITION  POINT 


R1  VERSUS  R2  AT  TRANSITION  POINT 


R1  VERSUS  R2  AT  TRANSITION  POINT 


11=23/13*6  =20HM/!0Ha‘(03a)0f= 


R1  VERSUS  R 2  AT  TRANSITION  POINT 


I 


Figure  E.7  Rj  vs  Rt  at  the  first  trans.  points,  for  P«  6°,  p|  p->  =  0  90,  Cj  c 


80 


RH01/RH02=.8,C1/C2=1.2 


=20HH/!0HH 
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